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This rapporteur review covers selected results presented in the Parallel Session HEA2 
(High Energy Astrophysics 2) of the 10th Marcel Grossmann Meeting on General Rel- 
ativity, held in Rio de Janeiro, Brazil, July 2003. The subtopics are: ultra high energy 
cosmic ray anisotropies, the possible connection of these energetic particles with powerful 
gamma ray bursts, and new exciting scenarios with a strong neutrino-nucleon interaction 
in the atmosphere. 



1. Introduction 

Since the early 60 's several ground-based experiments have observed extensive air 
showers, presumably triggered by ultra high energy cosmic rays (UHECRs) inter- 
acting in the upper atmosphere.l^ The highest primary energy measured thus far 
\s E ^ lO^'^'^ eV,^ corresponding to a center-of-mass energy -y/s = ^j2mpE ^ 
750 TeV, where is the proton mass. The interest in the origin of these particles 
is twofold: there is not only the intellectual curiosity about unknown properties of 
powerful astrophysical scenarios, but also the possibility to probe new physics at 
energies beyond the reach of any foreseeable man-made experiments. 

Theoretically, one expects the CR spectrum to fall off somewhat above 10^" eV, 
because the particle's energy gets degraded through interactions with the cosmic mi- 
crowave background (CMB), a phenomenon known as the Greiscn-Zatsepin-Kuzmin 
(GZK) cutoff.'"' Unfortunately, as one can see in Fig. ^ the most recent measure- 
ments by the HiReJ^l and AGAS^ experiments are apparently in conflict, if only 
statistical errors are taken into account, and the source of the difference remains 



'Rapporteur review of HEA2 - 10th Marcel Grossmann Meeting on General Relativity — 
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Figure 1. Data on the upper end of the cosmic ray energy spectrum with their statistical error 
bars, (HiResPAGASAP Fly's EyeP Haverah Parkland SUGArIHI) 

unknown. However, if one takes the systematic uncertainties in the energy measure- 
ments into account, one finds that both data sets are mutually compatible on the 
2 a level!^ Attempts to explain the AGASA data with a homogeneuos population of 
astrophysical sourc es that ini ects power-law distributions of C Rs g ive unacceptable 

(see, e.g., Refs. j^ l I j). On the other hand, an analysis - of the combined 
data reported by the HiRes, the Fly's Eye, and the Yakutsk collaborations is sup- 
portive of the existence of the GZK cutoff at the > Scr (> 3.7(7, depending upon 
the extrapolated energy spectrum) level.'' The deviation from GZK depends on the 
set of data used as a basis for power law extrapolation from lower energies. One 
caveat is a recent clainJi^ that there may be technical problems with the Yakutsk 
data collection. In view of the low statistics at the end of the spectrum and the 
wide variety of uncertainties in these experiments, perhaps the rational thing to do 
is to wait for more data, conservatively arguing that the jury is still out. 

In this Parallel Session we saw many thorough reviews covering all the most 



''This evidence disappears, however, if one assumes that UHECRs are protons and excludes nearby 
( < 50 Mpc) sources from the otherwise homogeneous distributionli^In this case even the HiRes-1 
data are incompatible with the GZK cutoff on the 3 o- levelESnU 
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interesting and timely topics in CR physics. In this rapporteur summary we cannot 
do justice to aU the presentations.*^ Priority will be given to two intriguing scenarios, 
which pose possible explanations of the data. 

2. Anisotropies in UHECRs 

At the highest energies, the arrival directions of CRs are expected to begin to reveal 
their origins. If the CR intensity were isotropic, then one should expect a time- 
independent flux from each direction in local detector coordinates, i.e., declination 
and hour angle. In that shower detected with local coordinates could have 

arrived with equal probability at any other time of a shower detection. For any 
point of the celestial sphere, the expected shower density can be estimated if the 
exposure in each direction can be obtained. This implies that celestial anisotropies 
can be easily discerned by comparing the observed and expected event frequencies 
at each region. 

For experiments with 100 % duty cycle, continuous operation in solar time for 
several years leads to a uniform observation in right ascension. Therefore, one of the 
conventional methods to search for any global anisotropy is to apply the Linsley' 

jni 

harmonic analysis to the full sky cosmic ray distribution, i.e., determine the ampli- 
tude and phase of the m**^ harmonic by fitting the right ascension distribution of 
events to a sine wave with period 27r/m. 

There is a remarkable agreement among several experiments favoring a signif- 
icant anisotropy (encoded in the first harmonic amplitude) around 10^* eV from 
the general direction of the Galactic Plane (GP). Specifically, the AGASA experi- 
ment has revealed a correlation between the arrival direction of CRs (with energy 
~ IQis eV) and the GP at the 4 a levell^ The GP excess, which is roughly 4% of 
the diffuse flux, is mostly concentrated in the direction of the Cygnus region, with 
a second spot towards the Galactic Center (GC)."'^^ Evidence at the 3.2 cr level for 
GP enhancement in a similar energy range has also been reported by the Fly's Eye 
Collaboration!2ii Interestingly, the full Fly's Eye data include a directional signal 
from the Cygnus region which was somew hat lost in an unsuccessful attempt to 
relate it to 7-ray emission from Cygnus Finally, the existence of a point-like 

excess in the direction of the GC has been confirmed via independent analysis of 
data collected with the SUGAR experiment 

For the ultra high energy (> 10^^'^ eV) regime, all experiments to date have 
reported no departure from isotropy in the first harmonic amplitudeP^ This does 
not imply an isotropic distribution, but it merely means that available data are too 



''A scen ario in wliicli UHECRs are able to break tlie GZK barrier was presented by She-SlienH 
XueEni 

"^For the Fly's Eye data-sample the first harmonic amplitude is computed using weighted showers, 
because it has had a nonuniform exposure in sideral time. A shower's weight depends on the hour 
of its sideral arrival time, and the 24 different weights are such that every time bin has the same 
weighted number of showers. 
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sparse to claim a statistically significant measurement of anisotropy. In other wor ds, 
there may exist anisotropics at a level too low to discern given existing statistics!^ 
The right ascension harmonic analyses are completely blind to intensity vari- 
ations which depend only on declination. Combining anisotropy searches in right 
ascension over a range of declinations could dilute the results, since significant but 
out of phase "Rayleigh vectors" from different declination bands can cancel each 
other out. Moreover, the analysis methods that consider distributions in one celes- 
tial coordi nate , while integrating away the second, have proved to be potentially 
misleading!^ An unambiguous interpretation of anisotropy data requires two in- 
gredients: expo sure to the full celestial sphere and analysis in terms of both celestial 
coordinates!^^ 

The first full sky search for large scale anisotropics in the distribution of arrival 
directions of CRs with energy > lO"'^^'^ eV was reported in this Parallel Session by 
John Swain. ' Data from the SUGAR and AGASA experiments, taken during a 
10 yr period with nearly uniform exposure to the entire sky, show no departures 
from either homogeneity nor isotropy on angular scale greater than 10°. 

In this full-sky anisotropy search, the intensity distribution of the set of = 99 
arrival directions 

1 ^ 1 

was conviniently expanded in spherical harmonics (Ygm) 

oo I 

J(n)=^ ^ ae„iYi,n{n) , (2) 

going into the multipole expansion out to £ — 5. Here, luj is the relative exposure 
at arrival direction Uj and A/" is the sum of the weights ojj^. The coordinate 
independent total power spectrum of fluctuations, 

cii) - t «L , (3) 

m——i 

is consistent with that expected from a random distribution for all (analyzed) mul- 
tipoles, though there is a small {2a) excess in the data for £ = sP^To give a visual 
impression of the level of homogeneity and isotropy in existing data, in Fig. |21 we 
show the intensity distribution as seen by AGASA and SUGAR experiments. 



3. UHECRs from GRBs 

In this section, arguments for the origin of UHECRs from gamma ray bursts (GRBs) 
are reviewed. This line of enquiry has led to a complete model for CRs originating 
from supernovae (SNe) and GRBs in our Galaxy and throughout the universe ;^^ 
which is summarized here. 
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Figure 2. UHECR intensity in arbitrary units (equatorial coordinates) as seen by the AG AS A 
and SUGAR experiments. 

The connection bet weeen G RBs and UHECRs was first made on the basis of 
an intriguing coincidencJ^OIS beween the power required to sustain the measured 
flux of super-GZK ( ^ 10^° eV) CRs against photohadronic energy losses and the 
local time- and space-averaged hard X-rays/soft 7-ray luminosity of GRBs. This 
luminosity density is estimated to be w 10'*'* ergs Mpc""^ yr~*- Thus GRBs have, in 
principle, sufficient energy to power the UHECRs and post-GZK CRs. Moreover, 
many GRB sources are found within the GZK radius, and CRs with energies ^ 10^° 
eV can be accelerated by the relativistic shocks formed in GRB explosions .E2 The 
hypothesis of a GRB/UHECR association points to a closer connection between 
SNe and CRs that could provide a complete solution to the problem of CR origin. 



3.1. CRs from Supernovae 

Even though the controversy surrounding the origin of the UHECRs has generated 
much interest, it should be noted that the much older problem of the origin of the 
CRs is itself not solved. Cosmic rays with energies from GeV/nucleon up to hun- 
dreds of TeV are widely thought to be accelerated by supernova remnant (SNR) 
shocks. Yet the prediction that SNRs should be luminous 7-ray sources and dis- 
play the characteristic 70 MeV vr'^ decay emission feature from hadronic interactions 
was not confirmed by the EGRET instrument on the Compton Observatory. Nev- 
ertheless, there is statistical evidence that SNRs are associated with unidentified 
7-ray sources.^I There is also clear evidence for a 7r° decay feature in the diffuse 
galactic 7-ray background, even if the spectrum is harder than would be expec ted 
if CRs throughout the Galaxy have the same spectrum as those observed locally.l^ 
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There is also some e vid ence from Cas A that the TeV emission originates from 
hadronic interactions^! though at a much lower level than would be expected if 
SNRs accelerated CRs with ^ 10 % efficiency. 

These clues suggest that the solution to the CR problem below the knee of the 
CR spectrum at « 3 PeV and, indeed, the problem of the origin of CRs at all 
energies, might be a consequence of the diversity in the types of supernovae (SNe). 
SNe are hardly uniform: the simplest separation is between the explosively burning 
white dwarfs (Type la), and the core co llap se SNe (Type II and Types lb and Ic). 
Moreover, the ejecta speeds vary greatlypfi' from a few thousand km per second in 
Type II SNe with massive H envelopes, to speeds of « 10^ km/s in Type Ib/c SNe 
where the H and He envelopes have been lost either to stellar winds or to Roche 
lobe overflow to a binary companion. 

The discovery that the GRBs which are observed today took place at cosmolog- 
ical distances led to the development of the relativistic fireball/blast- wave mode. 
that generalizes the theory of SNe to stellar explosions with relativistic ejecta. The 
coasting Lorentz factors of GRB outflows reach values of 10^ - 10^. Associated 
with this discovery are other far-reaching observational result J2SI that tie GRBs to 
a subset of SNe: GRBs are associated with massive stars in star-forming galaxies; 
optically dark GRBs could be due to extreme reddening from large quantities of 
dust and gas, as found in molecular cloud complexes where massive stars are born; 
and delayed reddened excesses in the late-time optical afterglow light curves appear 
to be SN emissions, and can often be fit with a template light curve of the Type 
Ib/c SN 1998bw. 

A final important discovery is that GRB outflows are in the form of highly col- 
limated jets. The evidence for beaming is inferred from beaming breaks in optical 
afterglow light curves that occur when the Doppler cone of the decelerating rela- 
tivistic outflow is about the same size as the jet cone. From the beaming breaks, 
one can infer that the mean solid angle subtended by a GRB jet is about l/500th of 
thefufl sky.EniAs a consequence, there are many hundreds of GRBs events that take 
place for every one that is observed. By performing the statistics of GRB sources, 
one finds that the rate of GRBs in the Milky Way reaches ^ 10% of the rate of 
Type Ib/c SNe. All these lines of evidence suggest that GRBs are a species of SNe. 

3.2. CRs from GRBs 

What does this mean for a complete model of cosmic rays? It is not sufficient to pro- 
vide an explanation for CRs with energies ^ 10^* eV without at least speculating 
on the origin of CRs with energies from the knee to the ankle. Many explanations 
have been suggested, including CR production from pulsars or extremely energetic 
SNe in the Galaxy, acceleration at the Galactic wind termination shock, or extra- 
galactic models where CRs above the knee diffuse into the GalaxyP^ A complete 
explanation based on acceleration in SNe shocks would tie our understanding of 
CR acceleration by SNe at GeV ~ TeV energies with an origin of UHECRs in the 
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relativistic shocks associated with GRB SNe. 

In the model of Ref. P|, high-energy cosmic rays (HECRs), defined as those 
with energies ^ 10^^ eV, are assumed to be accelerated at the shocks produced by 
the Type Ib/c SNe that collapse to form GRBs. From the evidence for beaming and 
the association of GRBs with star-forming galaxies like the Milky Way, GRB events 
are estimated to occur once every 3000 - 10000 yrs in the Galaxy. Relativistic shock 
acceleration in GRB blast waves are assumed to inject power-law spectra of ions 
from a minimum energy E'min ~ 10^"' eV to a high-energy cutoff i?max ^ 10^° eV. 
HECRs injected in the Milky Way diffuse and escape from our Galaxy, and UHE- 
CRs with energies ^ 10^^ - 10^* eV that have Larmor radii comparable to the 
size scale of the halo escape directly from the Milky Way and propagate almost 
rectilinearly through extragalactic space. By the same token, UHECRs produced 
from other galaxies can enter the Milky Way to be detected. UHECRs formed in 
GRBs throughout the universe travel over cosmological distances and have their 
spectrum modified by energy losses, so an observer in the Milky Way will measure 
a superposition of UHECRs from extragalactic GRBs and HECRs produced in our 
Galaxy. 

The model of Ref. fits the measured KASCADePH spectra of HECRs in the 
knee region, and the HiRes data^ at ultra-high energies. The fits imply the HECR 
injection spectral index and luminosity function of GRBs. The results show that 
GRBs must be strongly baryon-loaded, with the testable prediction that GRBs will 
produce a detectable number of high energy neutrino showers in a km-scale neutrino 
detector such as IceCube. 

3.3. GRB Model for HE and UHECRs 

High energy CRs from a GRB in the Galaxy are assumed to propagate diffusively as 
a result of pitch-angle scattering with magn eto- hy drodynamical (MHD) turbulence 
superposed on the Galactic magnetic field.l^E3 xhe Larmor radius tl of a CR 
propagating in a magnetic field of strength = l/^G is = A'^e/{ZB^G) pc, 

where 7 = 10^76 is the Lorentz factor of a CR with atomic mass A and charge 
Ze. The mean- free-path A between pitch-angle scatterings of a CR with Larmor 
radius is assumed to be inversely proportional to the energy density in the 
MHD spectrum at wave-number k ^ r^^- A two-component turbulence spectrum 
is assumed with wave- number index q = 3/2 for a Kraichnan-type spectrum at 
small wave number, and index q — 5/3 for a Kolmogorov-type spectrum at large 
wave numbers. The two components give a Z-dependent break in the scattering 
mean- free-path A at energy ^^^(PeV) = ZB^Qb^ci where we take B^q = 3, and find 
that 6pc = 1.6 is the wavelength in parsecs of the MHD waves where the spectrum 
changes from Kraichnan to Kolmogorov turbulence. The injection of turbulence 
at large wave numbers is probably due to SN explosions in the Galaxy; the small 
wave number turbulence could be a consequence of halo/disk interactions (e.g., 
high- velocity clouds passing through the Galaxy's disk). 
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Figure 3. Solid curve shows the value of (In A) calculated from the model of Wick, Dermer, and 
Atoyan (2004), where a single GRB in the Galaxy releases lO^'^ ergs in the form of a power-law 
spectrum of CRs with energies between 10^* eV and lO'^'^ eV. In the model, a Galactic GRB 
occurred 2.1 X 10^ years ago at a distance of 500 pc from Earth. The CRs isotropically diffuse 
via pitch-angle scattering with an energy-dependent mean-free-path A determined by the MHD 
turbulence spectrum. The dashed curves outline the expe cted range of uncertainty in (In A) from 
the experimental data and our knowledge of particle physicM!il The effect on (In A) of extragalactic 
CRs at energies 10^^ is not included. 



The diffusion radius rdif = 2 -^Z Act/3. Wiien r ^ rdif, the HECR difFerentiai 
number density n{'-f;r,t) oc t~^/^ x 'y^P^3(2-g)/2 rpj^^ measured spectrum is steep- 
ened by I (2 — q) units because the diffusion coefficient oc A oc 7^^^ for an 
impulsive source.ESI An injection spectrum with p = 2.2, as expected from relativis- 
tic shoclc acceleration, gives a measured spectrum nz,Ai"f',r,t) oc 7~*, with s = 2.7 
at <C Ez and s — 2.95 at E ^ Ez- These indices are similar to the measured 
CR indices below and above the knee energy. 

FigurcOshows the predicted mean logarithmic mass (In A) for CRs which diffuse 
from a GRB that occurred 2.1 x 10^ years ago at a distance of 500 pc from Earth. 
For comparison, the range of (In A) derived from direct measurements using a phe- 
nomenological model of particle interactions is shownP^ As can be seen, the model 
gives a good fit to the data at energies <C 10"'^^ eV, which is expected because the 
model also provides good fits to the individual spectra of CR ions measured with 
KASCADE. The model displays an excess abundance of heavy ions at higher ener- 
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gies, but the calculated value of (In A) shown here does not include contributions 
from extragalactic CRs, which become important above ^ 10^^ eV (see Fig. QJ. 
If the extragalactic component is primarily CR protons, then this discrepancy is 
resolved. 

Good fits to knee data were obtained with a break in the turbulence spectrum 
at wave numbers corresponding to « 1 pc. Composition enhancements by a fac- 
tor of 50 and 20 for C and Fe, respectively, over Solar photospheric abundances 
are also required. The likelihood for such an event is reasonable, and the corre- 
sponding anisotropy of the CRs from a single source is shown to be consistent with 
observations El 

UHECRs produced by GRBs throughout the universe lose energy adiabatically 
during the expansion of the universe, and through photo-pair and photo-pion pro- 
duction on the CMB. The loss processes produce features in the UHECR flux 
from distant {z ^ 1) GRB sources at characteristic energies ~ 4 x IQ-'^* eV and 
5 X 10^^ eV from photo-pair and photo-pion processes, respectively. Comparison 
of the integrated and evolved UHECR spectrum with data depends on how the 
luminosity density of GRB sources evolve through cosmic time. In view of the evi- 
dence that GRBs are associated with massive stars, the GRB luminosity density is 
assumed to be proportional to the star formation rate (SFR) history of the universe 
as traced by the blue and UV luminosity density of distant galaxies. The lowest 
possible SFR is expected to follow the evolution of the observed optical/UV lumi- 
nosity density, and a higher rate to follow this rate with allowance for extinction 
corrections.'^ 

Figure 0] shows the best fit model to the all-particle spectrum from below the 
knee of the CR spectrum to the highest energies. The best fit model to the data has 
an injection index p = 2.2, a high-energy exponential cutoff energy of 10^° eV, and 
a SFR which follows the higher extinction-corrected rate. The transition between 
galactic and extragalactic CRs is found in the vicinity of the second knee (10^^-^ eV), 
consistent with a heavy-to-light composition change .ESI The ankle (10^*-^ eV) is in- 
terpreted as a suppression from photo-pair losses, analogous to the GZK suppression 
at higher energies. 

Fits to the data define the local luminosity density of GRBs for produc- 
ing HECRs. The GRB HECR luminosity density is found to be 9^ 70 x 
10''''' erg Mpc~'^yr^^, so that this model implies that GRBs inject considerably 
more energy in the form of nonthermal hadrons than electrons. Thus GRB blast 
waves are baryon-loaded by a factor « 60 - 100 compared to the energy radiated by 
GRBs in the form of hard X-ray and soft 7-ray emission. For the large baryon load 
required for this model, calculations show that 100 TeV - 100 PeV neutrinos could 
be detected several time s per y ear from all GRBs with kilometer-scale neutrino de- 
tectors such as 

IceCube.™! 

Detection of even 1 or 2 neutrinos from GRBs with 
IceCube or a northern hemisphere neutrino detector would unambiguously demon- 
strate the high nonthermal baryon load in GRBs, and would provide compelling 
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Figure 4. Fit of the model of Ref. j^^j to the KASCADE (< 2 X 10^^ eV), HiRes-II Monocular 
(circles), and HiRes-I Monocular (inverted triangles) data, assuming that GRBs inject a power 
law distribution of particles with index p = 2.2 from lO'^* eV to a maximum energy with an 
exponential cutoff at -Emax = lO'^" eV. A minimum routine was used to find the best fit model 
which has a Xr — 1-03. This fit implies that the transition from Galactic to extragalactic CRs 
occurs near the second knee at 10^^'® eV, and that the ankle (lO'^®'^ eV) is a feature associated 
with photo-pair production. The heavy dark points shows the best fit model, and the galactic and 
extragalactic components are shown separately. 



support for this scenario for the origin of CRs. 

Finally, we note that these calculations have implicitly assumed that the Fly's 
Eye/HiRes measurements of the UHECR spectrum are correct. If a post-GZK 
excess is confirmed, then physics beyond the standard model could be required, as 
discussed in Sec. ^ 



3.4. The Last GRB in the Galaxy and the CR Excess from the GC 
As discussed in Sec. the AGAS^Pl and the SUGAbISI 

experiments observed 

an excess of CRs from the direction of the GC. The SUGAR data suggest a point 
source to within their spatial resolution, while AGASA shows an extended source. 
The excess starts to be significant around 10^^'^ eV, peaks near 10^^ eV, and cuts off 
sharply at about 10^^'^ eV. The flux of the excess particles represents a luminosity 
of particles beyond 10^^ eV of about 4 x 10'^'^ erg/s. The GC anisotropy is then 
very suggestive of neutrons as candidate primaries, because the directional signal 
requires relatively-stable neutral primaries, and time-dilated neutrons can reach the 
Earth from typical Galactic distances when the neutron energy exceeds 10^* eV. 

A novel hypothesis to explain the GC excess involving GRBs was presented in 
this Parallel Session by Gustavo Medina Tanco.E^ln order to estimate the remaining 
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traces of any CR activity produced by the last GRB in the galaxy, one has to take 

into account several considerations: 

(i) The UHECRs escaping the GRB fireball 



are mostly neutrons, because protons are captive in the magnetic field and suffer 
extensive adiabatic losses on the way outl^^P Some of these neutrons will decay 
into protons within the GC thin disk-like (r ~ 3 kpc) region of high interstellar 
medium density and high star formation rate. The population of secondary protons 
would then be captured by the strong i?- field near the GC, attaining diffusion with 
a residence time scale of about T ^ 10^ yr. At the end of this time, about 1/300 
protons are able to avoid leakeage. The trapped protons, 



can then be turned back into neutrons by interaction with nuclei in the interstellar 
medium with probability of 5 x 10"^. Here, m„ and t are the neutron mass and 
lifetime, respectively. 

(ii) The formation of the n ^ p reservoir depends on the GRB rate in the inner 
Galaxy times the probability that a GRB jet points more o less along the direction 
of the GP. The latter is estimated to be about 50 %. 
(in) The total CR production by a single GRB is w 10^^ erg.E^ 

Putting all this together, the observed anisotropy in the direction of the GG 
can be easily fitted by the neu trons produced in the GRB reservoir, that ultimately 
travel unscathed to Earth.l^ Arguably, if the anisotropy messengers are neutrons, 
then those with energies below 10^^ eV will decay in flight, providing a flux of cosmic 
antineutrinos above 1 TeV that should be observable at kilometer-scale neutrino 
telescopes.^ 

3.5. UHECRs, EMBHs, and GRBs 

An alt erna tive explanation for the GRB phenomenon has been recently put 
forward.E3 In this model the GRB explosion is due to the rapid discharge of an 
electromagnetic black hole (EMBH) . This chain reaction occurs when a massive star 
collapses into an EMBH able to produce a dyadosphere, i.e., the EMBH reaches the 
critical field for the Schwinger process, £c = mlc^/he ~ 1.3 x 10^^ V/m.l^ The 
pairs created through vacuum polarization then promptly annihilate into the GRB. 
This model explains the time variability, the spectra, and the GRB afterglows to a 
very high level of accuracy.l^ Based on the EMBH-GRB hypothesis, in this Parallel 
Session Alvise Mattel^ presented a mechanism to accelerate the ionized hydrogen 
atoms surrounding the death star to ultra high energies. 



No{E > 10^^ eV) - 10-2 No{E > E^^^) 



(4) 




(5) 



"^We remind the reader that the differential injection spectrum of GRBs oc E 
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Figure 5. Limits on charge and mass ratio for UHECR emission. EMBHs lying in the shaded zone 
develop a dyadosphere. EMBHs above the dotted line cannot activate the ionization mechanism. 
EMBHs on the right hand side of the dashed-dotted lines produce particles at lower energies 
between 10^^ eV and lO'^^ eV. The vertical solid lines indicate the maximum energy limit for proton 
acceleration, A, B, and C corresponds to Smax > 10^^ eV, Emax > 10^'^ eV, and Emax > 10^^ eV; 
respectively. (This figure is courtesy of Alvise Mattei). 

In this set-up protons would be accelerated outward in a very short time along a 
straight line by the induced electric field. Larmor losses would overcome gains if the 
energy increases as dE/dx ~ 10^^ MeV/m. Thus, protons do not suffer catastrophic 
energy losses if the EMBH cannot produce a dyadosphere. In terms of the mass 
AfsH and the charge Qbh of the BH, the condition for dyadosphere formation near 
the BH horizon (of radius rn), ^(t'h) < ^c, is found to be 

/ , \ -2 

< 6 X 10^^ 



1 + 



(6) 



where /i = M-byi/Mq and ^ = QBH/Omax- In order to produce a baryon reservoir, 
the electromagnetic field has to be able to ionize hydrogen atoms, i.e., a ionization 
potential of 13.6 V must be active in a distance of about a Bohr radius. This 
condition can be re-written as f (j'h) ^ 2.5 x 10^^ V/m, or equivalently 

^ < 2.9 X 10^^ e f 1 + ^/I-SA ^ ■ (7) 



The limits impossed by these constraints are summarized in Fig. 

In summary, an EMBH would lose its charge by emitting UHECRs if the BH 
does not reach the the critical charge-to-mass ratio given in Eq. This type of 
BH has been classified as pure or U-EMBH. The BH population that has recently 
exploded as a GRB ends up near the dyadosphere zone with residual charge, and 
so UHECR emission is still possible until complete discharge. 
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Akeno + AGASA TtI 




E [eV] 

Figur e 6. Ultraliigh energy cosmic ray data witli their statisti cal errors f to o: c ombination of 
AkenJMl and AGAS^El data; bottom: combination of Fly's EyfE^ and HiReJHHI data) and the 
predictions arising from a power-law emissivity distribution jSJ corresponding to sources which 
are uniformly distributed at cosmological distances. The best fits between E— = lO^^ '^ eV and 
= lO^'^ eV are given by the solid lines and correspond to the indicated values of the parameters 
a and n in the source emissivity distribution. The 2-sigma variations corresponding to the minimal 
(dotted) and maximal (dashed) fluxes are also shown. Other parameters of the analysis were 
E^ax = 3 X 10^1 eV, Zniin = 0.012, and 2max = 2. From Ref. Plj. 



4. Strongly Interacting Neutrinos 

We have noted in the Introduction that if the highest energy cosmic rays are nu- 
cleons (or nuclei), if their sources are indeed uniformly distributed at cosmological 
distances, and if their injection spectra are power-laws in energy - a reasonable 
assumption, in view of the measured spectrum in Fig. 0] which appears to be ap- 
proximately of (broken) power-law type over many order of magnitude in energy - 
then their total flux arriving at Earth should show a pronounced drop above the 
GZK cutoff Egzk w 4 X lO^^ eV (1 x 10^° eV, for nuclei). This is due to the fact 
that, above this energy, the universe becomes opaque to high energy nucleons (and 
nuclei), due to inelastic hadronic scattering processes with the CMB photons. The 
GZK cutoff is, however, not seen in the data, at least not in a significant manner 
(cf. Fig. El . Correspondingly, the events above 10^" eV in Fig. IHl should originate 
from small distances below 50 Mpc, the typical interaction length of nucleons above 
Eqzk- However, no source within a distance of 50 Mpc is known in the arrival 
directions of the post-GZK events. The basic puzzle is: if there are no large inter- 
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veiling magnetic fields and the sources of ultraliigli energy cosmic rays are indeed 
at cosmological distances, liow could they reach us with energies above 10^° eV? 

At the relevant energies, among the known particles only neutrinos can prop- 
agate without significant energy loss from cosmological distances to us. This fact 
leads naturally to scenarios invoking hypothetical - beyond the Standard Model 
of elementary physics - strong interactions of ultrahigh energy cosmic neutrinos,!^ 
whose modern incarnation we review in this section. '^^ 

Such scenarios are based on the observation that the flux of neutrinos orig- 
inating from the decay of the pions produced during the propagation of nucle- 
ons through the CMbSS.'SI] _ the cosmogenic neutrinos - shows a nice agreement 
with the observed UHECR flux above Eqzk- Assuming a large enough neutrino- 
nucleon cross-section at these high energies, these neutrinos could initiate extensive 
air showers high up in t he atmosphere, like hadrons, and explain the existence of 
the post-GZK events.Ei This large cross-section is usually ensured by new types of 
TeV-scale interactions beyond the Standard Model, such as arising through gluonic 
bound state leptons,IS2 through TeV-scale grand unification with leptoquarks,E3 
through Kaluza-Klein modes from compactified extra dimensions,!^ (see, however, 
Ref. [ J), or throughp-brane production in models with warped e xtra dimensions 
(see, however, Ref. p^); for earlier and further proposals, see Ref. and Ref. 
respectively. 



In Refs. [EIEn|^ a detailed statistical analysis of the agreement between obser- 
vations and predictions from such scenarios was presented. Moreover, an example 
was emphasized which - in contrast to previous proposals ~ is based entirely on 
the Standard Model of particle physics. It exploits non-perturbative electroweak 
instanton-induced processes for the interaction of cosmogenic neutrinos with nucle- 
ons in the atmosphere, which may have a sizeable cross-section above a threshold 
energy Eth = 0{{^TTmw / awY) / {'^'mp) = 0(10^^) eV, whe re tom^ de notes the W- 
boson mass and aw the electroweak fine structure constant .EHHIEl 

The scenario is based on the assumption of a power-law emissivity distribution 
corresponding to uniformly distributed sources and thus quite consistent with the 
GRB model of CR origin from Sec. 18.31 The emissivity is defined as the number of 
protons per co-moving volume per unit of time and per unit of energy, injected into 
the CMB with energy Ei and characterized by a spectral index a and a redshift [z] 
evolution index n, 

Lp = Jo (1 + zf 0{E^,^ - E,) 0{z - Z,,in) 0(z,,ax " z) . (8) 

Here, jo is a normalization factor, which will be fixed by the observed flux. The 
parameters -Emax and Zmin/max have been introduced to take into account certain 
possibilities such as the existence of a maximal energy, which can be reached through 
astrophysical accelerating processes in a bottom-up scenario, and the absence of 
nearby/very early sources, respectively. They turn out to be quite insensitive to 
the specific choice for i^max, -Zmin, and Zmax7 within their anticipated values. The 
main sensitivity arises from the spectral parameters a and n, for which 1 and 2 a 



1,2008 23:19 WSPC/Trim Size: 9.75in x 6. 5in for Proceedings 



rio'crS 



15 




Figure 7. Ultrahigh energy cosmic ray data (Akeno + AGASA on the upper panel and Fly's Eye 
+ HiRes on the lower panel) and their best fits (solid) within the electroweak instanton scenario, 
for -£/niax — 3 X 10^2 eV, ^min = 0.012, 

^max — 2, consisting of a pro ton component (dotted) plus 
a cosmogenic neutrino-initiated component (dashed). From Ref. 



confidence regions have been determined. Note, that n ~ 3, ^max ~ 2 mimicks the 
SFR history of the universe mentioned also in connection with GRBs in Sec. 13.31 

After propagation through the CMB, the protons from Eq. (jSJ will have en- 
ergies below Eqzk, so they can well describe the low energy part of the UHECR 
spectrum. The cosmogenic neutrinos interact with the atmosphere and thus give 
a second component to the UHECR flux, which describes the high energy part of 
the spectrum. The relative normalization of the proton and neutrino fluxes is fixed 
in this scenario, so the low and high energy parts of the spectrum are explained 
simultaneousl y w ithout any extra normalization. Details of this analysis can be 
found in Ref. 

Figure □ shows the best fits for the AGASA and for the HiRes UHECR data. 
The best fit values are a = 2.68(2.68) and n = 2.65(2.9), for AGASA (HiRes), 
within the electroweak instanton scenario. One can see very nice agreement with 
the data within an energy range of nearly four orders of magnitude. The fits are 
insensitive to the value of E'max as far as one chooses a value above « 3 x 10^^ eV. 
The maximum injection ener gy however , is constrained by EGRET measurements 
of the diffuse gamma ray flux.L^ This is because the photons produced via 7r° decay 
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n 



Figure 8. Confidence regions in tlie a-n plane for fits to the Akeno + AGASA data (2-sigma (long 
dashed)) and to the Fly's Eye + HiRes data (1-sigma (solid); 2-sigma (short-dashed)), respectively, 
within the electroweak instanton scenario, for i?niax >3 X 10^1 eV, 
Ref. P|. 

are degraded in energy due to interactions with the universal radition backgrounds. 
It is noteworthy that the use of a 2-component model allows -Emax ~ 3 x 10^^ eV 
without violating the new EGRET boundsPi 

The shape of the curve between 10^^ eV and 10^^ eV is mainly determined by the 
redshift evolution index n. At these energies the universe is already transparent for 
protons created at z « 0, while protons from sources with larger redshift accumulate 
in this region. The more particles are created at large distances - i.e. the larger n is - 
the stronger this accumulation should be. In this context, on may note that the data 
seem to confirm the implicit assumption that the extragalactic uniform UHECR 
component begins to dominate over the galactic one already at « 10^^ eV. If one, 
alternatively, starts the fit only at 10^^'^ eV - corresponding to the assumption that 
the galactic component dominates up to this energy - one finds, however, also a 
very good fit, with a very mild dependence on n and the same best fit values for a, 
with a bit larger uncertainties. The peak around 4 x 10^^ eV in Fig. [3 shows the 
accumulation of particles due to the GZK effect. Neutrinos start to dominate over 
protons at around 10^" eV. 

Figure |H1 displays the confidence regions in the (a,n) plane for AGASA and 
HiRes. The scenario is consistent on the 2-sigma level with both experiments. 
For HiRes, the compatibility is even true on the 1-sigma level. It is important to 
note that both experiments favor the same values for a and n, demonstrating their 
mutual compatibility on the 2-sigma level (see also Ref. f^). If one ignores the 
energy uncertainty in the determination of the goodness of the fit, they turn out to 
be inconsistent. 
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Finally, let us emphasize that the same fit results are valid for all strongly in- 
teracting neutrino scenarios, as long as the neutrino-nucleon cross-section has a 
threshold-like behavior as in the case of electroweak instanton-induced processes, 
with a neutrino threshold energy < 4 x 10^^ eV and a cross-section > 1 mb above 
threshold. It is also important to note that the energy requirements on the sources 
of the primary protons are comparatively mild. To obtain a good fit, one needs 
-Emax > 3 X 10^^ eV. There are several astrophysical source candidates with CR- 
emission up to this energyP notably, as we discussed in Sec. 3, GRBs can provide 
the necessary conditions to accelerate protons to the required energies by conven- 
tional shock acceleration. Moreover, they naturally provide the power-law source 
emissivity distribution assumed in Eq. (jSjl. 

The predicted ultrahigh energy cosmic neutrino component can be experimen- 
tally tested at cosmic ray facilities by studying the zenith angle dependenc 1 1 ^5 | 
of the events in the range 10^* - 10^° eV and by analyzing possible correla- 
tions with distant astrophysical sources. Additionally, one can look for bumps 
in neutrino-initiated shower spectra at neutrino telescopes such as IceCube.l^ As 

laboratory tests, one may search for enhancements in (quasi-)elastic lepton-nucleon 
77 

scattering' or for signatures of QCD instanton-induced processes in deep-inelastic 
scattering, '^S' e.g. at HERA.E^ 

In summary, strongly interacting neutrino scenarios provide a viable and attrac- 
tive solution to the UHECR puzzle and may be subject to various crucial tests in 
the foreseeable future. 



5. Concluding Remarks 

We have reviewed some of the interesting ideas currently invoked to explain the 
origin of UHECRs. We expect future CR experiments, such as the Pierre Auger 
Observatory,*^ would clarify the confusing experimental situation and shed light 
on validity of models discussed in this report. 

The Pierre Auger Observatory is designed to measure the energy and arrival 
direction of UHECRs with unprecedented precision. It will consist of two sites, 
one in the Northern hemisphere and one in the Southern, each covering an area 
S « 3000 km'^ . The Southern site is currently under construction while the Northern 
site is pending. Once complete, these two sites together will provide the full sky 
coverage and well matched exposures which are crucial for anisotropy analyses. The 
base-line design of the detector includes a ground array consisting of 1600 water 
Cerenkov detectors overlooked by 4 fluorescence eyes. The angular and energy 
resolutions of the ground arrays are typically less than 1.5° and less than 20%, 
respectively. The detectors are designed to be fully efhcient out to ^max = 60° 
beyond 10^^ eV. In 10 yr of running the two arrays will collect « 4000 events above 
10^^'^ eV. Such statistics will enable us to solve the GZK puzzle. 
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